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Abstract. We study the causal structure of dynamical charged black holes, with a 
sufficient number of massless fields, using numerical simulations. Neglecting Hawking 
' radiation, the inner horizon is a null Cauchy horizon and a curvature singularity due to 

mass inflation. When we include Hawking radiation, the inner horizon becomes space- 
^ ' like and is separated from the Cauchy horizon, which is parallel to the out-going null 

<sJ '< direction. Since a charged black hole must eventually transit to a neutral black hole, we 

^ ] studied the neutralization of the black hole and observed that the inner horizon evolves 

into a space-like singularity, generating a Cauchy horizon which is parallel to the in- 
going null direction. Since the mass function is finite around the inner horizon, the 
inner horizon is regular and penetrable in a general relativistic sense. However, since 
the curvature functions become trans-Planckian, we cannot say more about the region 
', beyond the inner horizon, and it is natural to say that there is a "physical" space-like 

singularity. However, if we assume an exponentially large number of massless scalar 
fields, our results can be extended beyond the inner horizon. In this case, strong cosmic 
I censorship and black hole complementarity can be violated. 

OO 
O 

> ■ 

•'-^ . 1. Introduction 

X 

I The Reissner-Nordstrom black hole [1] has been a well-known solution since black hole 
physics began, see Figure 1. However, although the geometrical structure of the static 
charged black hole solution was well understood [2] , there still remained some interesting 
questions. 

For example, the Reissner-Nordstrom black hole has a time-like singularity, and this 
may imply that there exist causally undetermined regions. The boundary between the 
determined and undetermined regions is called the Cauchy horizon [3]. Some authors 
argued that there may be local effects like infinite blue shift [4] or semi-classical effects 
[5] along the Cauchy horizon and that these effects might imply that an observer will 
not be able to penetrate the Cauchy horizon. This would rescue the philosophical idea 
known as strong cosmic censorship, which states that no singularity is ever visible to 
any observer [3]. 
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In a pioneering work [6], Poisson and Israel discovered mass inflation as an 
important property of the inner Cauchy horizon. If there exists an energy flow along 
the inner horizon, and if an observer crosses it, he will feel infinite energy. This means 
that the inner horizon is unstable and now we understand that a realistic inner horizon 
must be interpreted as a kind of curvature singularity 

However, there still remains an important question: what will be the back-reaction 
of the instability on the causal structure? There has been some controversy on this 
issue. Some authors argued that the inner horizon^ is unstable and it may collapse to 
the singularity at the center [7]. Others said that, since the inner horizon is problematic, 
there should be no inner horizon in real situations due to pair-creation of charged 
particles [8] (see also [9]). Others argued that although the inner horizon has some 
problems, like a curvature singularity, it is stable and the inside structure seems to be 
safe, in the sense that the metric perturbation is finite and the tidal deformation is 
small enough [10] [11]. In other contexts, some authors argued that there are two kinds 
of singularities, strong and weak singularities [12], and that the inner horizon becomes 
a weak and null singularity [13], that is, an observer may hit the inner horizon first and 
then fall to the strong singularity deep inside of the black hole [14]. However, all of 
these ideas were based on a charged black hole without Hawking radiation. 

To figure out which opinion is correct and which is not, we need to understand the 
causal structure of charged black holes not only for static cases but also for dynamical 
cases which include the formation and the evaporation of the black hole. Of course, some 
previous authors, motivated by the information loss problem, considered the effect of 
Hawking radiation. They calculated the near extreme limit of a charged black hole 
using some approximations [15], or the Vaidya metric [16] [17], but these calculations 
could not include the mass infiation effect. The first successful implementation of mass 
infiation was in Hod and Piran [19] by using numerical simulations. 

Pioneering black hole simulations were done in [18] [19] [20] [21] [22] (see also [8]) whose 
methods we followed. In [19], the formation of a charged black hole from collapsing 
charged matter fields was studied, and a space-like singularity and a null inner Cauchy 
horizon were obtained (Figure 2). These results were confirmed and refined in [21] [22]. 
The first numerical simulation to include Hawking radiation for a charged black hole was 
done in [20], and it was observed that the outer and inner apparent horizons approach 
along the out-going null direction. 

In our simulation, we assume spherical symmetry and set up equations for 
Maxwell and scalar fields including Hawking radiation. After reproducing the results of 
previous authors, we did our own experiments about charged black holes with Hawking 
radiation and discharge. Especially, we focused on the causal structure and its physical 
implications which were not deeply discussed by previous authors. 

In this paper, we will answer the following questions: 

^ In this paper, we use the outer and inner horizons in a local sense, i.e. we will use "apparent" or 
"trapping" horizons [29]. 
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• What is the Penrose diagram for dynamical charged black holes with Hawking 
radiation and discharge? 

• Is the inner horizon stable or unstable? Is the inner horizon penetrable or not? 

• What are the implications for cosmic censorship and black hold complementarity? 

In Section 2, we describe our model for a dynamical charged black hole. In Section 3, we 
discuss the causal structure of dynamical charged black holes. In Section 4, we discuss 
the properties of the inner horizon. In Section 5, we discuss the implications for strong 
cosmic censorship and black hole complementarity. 

2. Model for a dynamical charged black hole 

We will briefly introduce the model used in this paper. For details, see Appendices 
Appendix A, Appendix B, and Appendix C. In this paper, we set G — c — kB — 47reo = 
1, and mpi = Ipi — Vh, where mp\ and lp\ are the Planck mass and length. 

We assume a complex massless scalar field (p which is coupled to the electromagnetic 
field [2]: 

>C = - + ieA„0) (0., - ieA,'^) - ^F^bF'"'', (1) 
where Fab = ^6;a — Aa;b, and e is the gauge coupling. We assume spherical symmetry 

ds^ = ^ o? (u, v) du dv + (-u, v) dVt^ , (2) 

where u is the retarded time, v is the advanced time, and 6 and are angular coordinates 
[23]. Spherical symmetry fixes = = 0, and we can choose the gauge Ay = giving 
= (a,0,0,0) [21]. 

To obtain the metric, we use the semi-classical Einstein equation: 

G,, = 8n (t^c, + (f « )) , (3) 

where the right hand side contains the energy-momentum tensor for classical collapsing 
fields and Hawking radiation (T^). For the expectation value of operators to 
have physical meaning on a classical background, the operator — (T^) of all 
possible quantum fluctuations should be sufficiently smaller than (T^), or roughly, the 
normalized dispersion [((T^)^) — (T^)^]/ (T^)^ should be sufficiently small. If we have a 
sufficiently large number N of massless fields, and if each field contributes independently 
to the energy- momentum tensor, then the normalized dispersion will decrease as 1/N. 
Therefore, using Equation (3) is justified in the large N limit. 

To include the Hawking radiation in the form of the renormalized energy- 
momentum tensor (Tj[^), we used 1 + 1 dimensional results [24] divided by 47rr^, 
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Figure 1. The Penrose diagrams of static charged black holes for Al > Q and M = Q. 
singularity 




Figure 2. Formation of a charged black hole with mass inflation, but without Hawking 
radiation. There is a null inner Cauchy horizon and a space-like singularity deep inside 
of the black hole. 



This is a reasonable model for the spherically symmetric case [20]. In Equation (4), a 
subscript of a function denotes a partial derivative, and 

P oc (5) 

where is the number of massless scalar fields generating the Hawking radiation and 
/pi is the Planck length. By changing P, we can tune the strength of quantum effects, 
and by changing at fixed P, we can tune the Planckian cutoff. 

For initial conditions which make the black hole, we take the charged field 
configuration at the initial surface to be 

0BH^^ y) = ^ sin2 ( vr^i^^ exp ( 27rz^^^ (6) 
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Figure 3. Contour diagram of r for (e = 0.1, P = 0, A = 0.25). This diagram confirms 
Figure 2. The outer horizon (ry = 0, red curve) grows in a space-hke direction, and 
the inner horizon (also, = 0) would be located at the v ^ oo limit. Here, spacing is 
1 for black contours and 0.1 for green contours. 



Figure 4. The contours of the mass function for log |m| = 0, 20, 40, 60, 80, 100, 200, 
300. 



for Vi < V < Vf and otherwise (f)^^{ui,v) = 0, where Ui = and f i = are the initial 
retarded and advanced time, and V{ = 20 is the end of the pulse in the initial surface. 
Finally, after fixing some parameters, we will be left with three free parameters: the 
gauge coupling e, the strength of Hawking radiation P, and the amplitude of the field 
A. We determine one specific simulation by choosing values for these three parameters 
(see Appendix B). 

3. The causal structure of dynamical charged black holes 

3.1. Formation of a black hole via collapsing matter fields 
The static charged black hole solution is well-known: 

ds^ = -[l + ^]de+[l + ^ dr^ + r^dn\ (7) 

From this metric we can draw a maximally extended causal structure giving the Penrose 
diagrams [2] in Figure 1. However, these diagrams are not true for dynamical situations. 
For example, the initial state should be generally fiat, and the final state will be fiat 
again in many cases. 

Then, what will happen if we collapse some charged matter in fiat space-time? 
One simple guess is the growth of outer and inner horizons, and the evolution of a 
time-like singularity [16]. However, because of mass inflation, which makes the inner 
horizon singular, it was guessed that there is a space-like singularity deep inside the 
black hole and a null inner horizon [10] [14] as in Figure 2. This idea was confirmed by 
some numerical simulations [19] [21]. Also, we reproduced the same result in Figure 3. 

To check whether the inner horizon becomes singular or not, we need to see the 
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behavior of the mass function. The mass function is defined by [26] 

(1 + ^+4!^), (8) 

where q{u,v) = 2r^at,/a^ is the charge (see Appendix A). Figure 4 shows the behavior 
of the mass function. It blows up exponentially as v ^ 00. This is a typical signature 
of mass inflation, m ~ exp(Kif), where Kj is the surface gravity of the inner horizon 
[6]. Then, some scalar quantities such as the Ricci scalar and the Kretschmann scalar 
(K = Ra/s-ysR""^^^ ~ rn^) become infinite as v ^ 00 [10]. Since we did not consider 
quantum effects, the Planckian cutoff of the scalar curvatures can be regarded as 00. 
Thus, the curvature singularity occurs at the v ^ 00 limit. 

Note that there is a transition region between the black curves and the green curves 
in Figure 3 where the decrease of the radial function begins to vary slowly. From 
Figure 4, we can see that mass inflation begins around there. This region corresponds 
to r ~ M - v'M^ - i.e. the inner horizon of the Reissner-Nordstrom metric, 
although it is not a real horizon. Its physical meaning and more detailed analysis will 
be included in a future work [27]. 

3.2. Evaporation and discharge 

The evolution of a charged black hole is determined by two quantum effects: Hawking 
radiation and charge pair creation. 

In a Reissner-Nordstrom metric, the rate of change of the mass depends on the 
Hawking temperature. From the Stefan- Boltzmann law, 

dt ^ " 2407r r% ' ^ ' 

where a — 7r^/60 is the Stefan- Boltzmann constant. 



r+ = M + ^M"^ - g2 (10) 
is the outer horizon, and 

_ 1 ^/M^^ 

^^"2^ ^^^^ 

is the Hawking temperature. 
The pair creation rate is 



e'E^ 



r ^ ^e-- (12) 



where 



Ec^ (13) 

e 

is the critical electric field, me is the electron mass, and e is the gauge coupling [31]. 
The electric field near the horizon is 

- (14) 
' + 
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Figure 5. Evolution of charged black holes, (a) If E+ ^ Ec, discharge dominates 
and the black hole approaches — E^- (b) If £'+ ^ Ec, discharge is exponentially 
suppressed. If Q/M <^ 1, Hawking radiation dominates and the mass will decrease 
via Hawking radiation, (c) If Q 3> -E^^, the black hole will approach and follow the 
extreme black hole (line Q = M). (d) If Q <C E~^ , the black hole will approach and 
follow E+ = Ec- (e) When the charge reduces to the last few quanta, Q ~ e, then 
Th ^ M^^ ^ TTie on the = Ec track, and the black hole will emit its final quanta 
of charge via Hawking radiation, (f) The final neutral black hole will lose mass via 
Hawking radiation. 



and the pair creation rate near the horizon is r+ = r{E^). Therefore, the rate of change 
of the charge becomes approximately 

dO o e2g2 

^ ~ r+V ~ r+r^ ~ (15) 

where V is the proper volume for pair creation [20]. 

Using these results, we can draw a schematic diagram for the evolution of a charged 
black hole as in Figure 5. 

(a) If E+ > Ec, then 

1 dM 1 (M2-Q2)2 g2g 1 



M dt M rl r+ Q dt' 

and hence pair-creation dominates Hawking radiation. Then the black hole will 
emit charge to approach the track E^ = E^. 



EM=^^{l + EcQ) {EcQ<l). (17) 

(b) If £■+ ^ Ec, discharge is exponentially suppressed and, if Q/M <^ 1, the black hole 
will lose mass via Hawking radiation. 

(c) If Q ^ E^^, then E^ ^ Ec still holds for the extreme black hole."*" Therefore, 
the black hole will approach an extreme black hole and the Hawking radiation also 

+ E^^ = e/TTTTig ^ lO'^^mpi ^ lO^Msun. To approach extremality, Q ^ lO^Afsun is required. 
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becomes suppressed. The black hole will then emit mass and charge maintaining 
extremality. 

(d) If Q ^ E^^, the black hole will approach the = track. The black hole will 
then start to emit charge and follow the track. 

(e) When the charge reduces to a few quanta, Q ~ e, then M ~ l/mg on the = E^ 
track, and the Hawking temperature becomes of the order of the electron mass, 
i.e. Th ~ ~ mg. Therefore, Hawking radiation will emit the final quanta of 
charge and the black hole will be neutralized [32]. 

(f) After the black hole is neutralized, it will lose mass via Hawking radiation. 
Thus, in terms of causal structures, we need to understand two important stages: 

• a non-extreme black hole evolves toward an extreme black hole via Hawking 
radiation, 

• the near extreme black hole evolves to a neutral black hole via discharge. 
We will observe these two stages independently. 

3.2.1. From a non- extreme black hole to an extreme black hole via Hawking radiation 
In this section, we investigate a charged black hole with Hawking radiation using the 
parameters (e = 0.1, P = 0.1, A = 0.25). Figure 6 and 7 show contour diagrams of 
r as well as the outer and inner horizons, which are defined by = 0. One can see 
that the inner horizon is space-like and approaches the out-going null direction as time 
V increases. The space-like inner horizon is different from Section 3.1, but agrees with 
[20]. 

In Figure 6, r increases beyond the inner horizon. This means that the r„ = 
contour is almost the same as the inner horizon r„ = 0. If we change the mass M and 
the charge Q, the existence and properties of the r„ = contour and the qualitative 
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Figure 7. Detailed plot near the outer horizon of Figure 6. This shows that the outer 
horizon (red curve) bends in a time-like direction. Here, spacing is 1 for black contours 
and 0.002 for green contours. 

Figure 8. r function along the green cutoff line of Figure 6, u ~ 18.18. The sharp 
point has r = \/P ~ 0.32, which is the central singularity. This confirms that the 
green cutoff line of Figure 6 is a Cauchy horizon. 




Figure 9. The Penrose diagram for Figure 6. 

behavior of the radial function can change. This behavior will be discussed in future 
work [34]. However, the r„ = horizon structures are qualitatively the same. 

One interesting feature is the green null cutoff line along u ~ 18.18 of Figure 6. 
We plotted the r function along this line in Figure 8, and it shows a sharp point at 
r = ~ 0.32, i.e. the central singularity. Although this cutoff line itself is regular, the 
region beyond cannot be calculated because of the singularity. This null section grows 
from the singularity, and we can interpret the null cutoff line as the Cauchy horizon by 
definition. Note that the inner horizon is separated from the Cauchy horizon. 
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Figure 10. Asymptotically calculated charge and mass during the neutralization. 

Figure 11. r function along the green cutoff line of Figure 12, v ~ 58.11. This 
confirms that the second green cutoff line is also a Cauchy horizon. 

When the mass of a black hole decrease via Hawking radiation, the inner horizon is 
space-like and the outer horizon is time-like [29]. Therefore, to connect the non-extreme 
Penrose diagram (Figure 2) to the extreme one (right of Figure 1), the only possible 
way is the time-like bending of the outer horizon, and the space-like bending of the 
inner horizon [15]. Therefore, this logical expectation is consistent with our numerical 
simulations. 

Finally, we can draw a schematic diagram in Figure 9.* 

3.2.2. From an extreme black hole to a neutral black hole via discharge In this section, 
we investigate a charged black hole with discharge. A transition from a charged black 
hole to a neutral black hole can be simulated by supplying opposite charged matter to 
the black hole: 

/idis / dis \ / dis \ 

for vf^^ < V < vf^^ and otherwise (j)'^^^{ui,v) = 0, where A'^''^ = 0.15, vf^^ = and 
ti^is = 120. Then the total scalar field initial conditions are 0i = (pf^ + As 
we slowly add this matter, the total charge decreases greatly, but the mass increases 
slightly, see Figure 10. We may regard this experiment as a good toy model for a general 
charged-neutral transition. 

In Figure 12, as the charge of the black hole decreases, there arises another null 
cutoff line in the u direction. As before, using Figure 11, one can see that this cutoff 
line is a Cauchy horizon. 

During the neutralization, the charged black hole evolves into a neutral black hole 
(left of Figure 13). Therefore, the region beyond the inner horizon should disappear, 

* Some previous authors pasted a near extreme black hole (left of Figure 1) to an extreme black hole 
(right of Figure 1) [15]. 
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Figure 12. Contour diagram of r. Note tliat, compared with. Figure 6, there is another 
Cauchy horizon along the u direction. 




Figure 13. The Penrose diagram for a transition from a near extreme charged black 
hole to a neutral black hole. 



and a space-like singularity should appear. In our simulation, the inner horizon evolved 
into the space-like singularity, and a Cauchy horizon formed. Finally, we can draw the 
causal diagram for this process (right of Figure 13). Our results are consistent with 
previous work [17]. 

3.3. Conclusion 

Now we can draw the causal structure of dynamical charged black holes by combining 
Figures 9 and 13 to obtain Figure 14. One can see the outer and inner horizons, as well 
as the Cauchy horizons. 
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4. Properties of inner horizons 

4.1. Penetrability of inner horizons 

The black hole was formed by the charged matter pulse 0^^, and to test penetrability, 
a second neutral pulse 0^°"^ was added to the initial conditions: 

/ipon / pen N 

0P-(ni, v) = ^ sin^ n ' (19) 

for v?^" < w < f^'^'^ and otherwise (f)^'^^{ui,v) = 0, where A^^^ = 0.01, ff"^"^ = 30 and 
■yP*^" = 40. This simulation is the same as Figure 6 between v = and f = 30, but 
modified from v = 30. The scalar field has initial conditions (pf^ + (pf'^'^. Note that 
the penetrating pulse has energy oc {A^^""^)"^ ~ 10~^ while the energy of the first pulse 
is oc ~ 0.06, and hence the perturbed energy is approximately 1/100 of the total 
mass. Since the equation for the scalar field. Equation A. 28, is linear in the field and 
its derivatives, and the energy of the penetrating pulse is sufficiently smaller than the 
energy of the initial pulse, we can trace the effect of (pf^ and 0?'^'^ independently, up to 
gravitational back reaction. 

Figure 15 shows the behavior of each pulse. In this diagram, one can see that 
the inner horizon has a barrier property; field configurations change drastically near 
the inner horizon. However, the inner horizon is penetrable; fields penetrate the inner 
horizon, and fiow beyond it. 
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Figure 15. Contour diagram of •\/47r3?((/)P*^"). It is affected by, but penetrates, the 
inner horizon. Here, spacing is 0.06457 for black contours and 0.002 for green contours. 

Figure 16. The contours of the Ricci scalar for log \R\ = 0, 20, 40, 60, 80, 100, 200, 
300, 400, 500. 

4-2. Mass inflation and the trans- Planckian problem 

Now we observe the curvature function for Figure 12. Figure 16 shows the contours 
of the Ricci scalar R. The Ricci scalar starts to blow up exponentially near the inner 
horizon, which is a typical signature of mass inflation [19] [21], but it does not diverge 
except at the central singularity r = y/P and it is still finite even on the Cauchy horizon. 

However, since we consider quantum effects, we should also consider the Planck 
scale. As the metric function a in Equation (2) decreases exponentially beyond the inner 
horizon, the mass function (~ 1/a^), the Ricci scalar (~ 1/a^), and the Kretschmann 
scalar (~ l/a^) blow up exponentially. When they become greater than the Planckian 
cutoff, the classical picture will break down. 

However, from Equation 5, the Planckian curvature cutoff, i?cutofr = lp\ oc N/P, 
is not defined until one has defined the number of massless fields A^. So, if we 
assume sufficiently large N, we can extend the classical picture to arbitrary large R. 
The required is of the order of the curvature R beyond the inner horizon, since 
R < RcntoS OC N/ P is required. The characteristic scale of R is of the order of the mass 
function M, which is of the order of the mass inflation factor exp K,i{u + v) [6], where Ki 
is the surface gravity of the inner horizon. K[ is of the order of 1/M and u and v are of 
the order of the lifetime and hence the required A^ is of the order of exp M^. 

If one can assume such a large A^, the whole region in Figure 14 becomes reliable 
except for the central singularity. The authors think that this would have important 
implications for cosmic censorship and black hole complementarity, see Section 5. 
However, if A^ is bounded in principle or such a large A^ spoils some assumptions 
of semiclassical gravity, we cannot avoid the trans-Planckian curvature problem. In 
this case, we should regard the trans-Plankian inner horizon as a kind of singularity 
(Figure 17).^ 

tt This is quite similar to some classical dilaton black hole models [33] , in which the inner horizon is a 
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Figure 17. The causal structure of dynamical charged black holes for a small number 
of massless fields. It is natural to think that there is a kind of physically singular 
region around the inner horizon. The region (red dashed curve) is space-like and has 
non-zero radius, while the central singularity has zero radius. 



5. Cosmic censorship and black hole complementarity 

Some authors have suggested that weak cosmic censorship, which states that all 
singularities are hidden within black holes, can be violated in charged black holes [16]. 
However, those authors used the Vaidya metric, in which the mass and charge functions 
are put in by hand, which can be unphysical. Weak cosmic censorship was not violated 
in our simulations (see also [21]). Whether strong cosmic censorship holds is less clear. 
Without Hawking radiation, the mass function diverges at the inner Cauchy horizon, 
and this prevents the violation of strong cosmic censorship. But if we include Hawking 
radiation, the mass function is finite at the Cauchy horizon, and if we assume sufficiently 
large N, all curvatures can be less than the Planckian cutoff, and hence strong cosmic 
censorship is violated. 

Finally, using our results, we investigate black hole complementarity [30] [35]. 
According to black hole complementarity, after the information retention time (when 
the area of the black hole decreases to half its initial value), an observer outside a 

physical singularity with non-zero radius, and the black hole approaches a stable extreme black hole. 
However, there are some differences: as discussed above, the inner horizon in our model is regular in a 
general relativistic sense, in that it is geodesically complete and penetrable. 
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Figure 18. A schematic diagram for the duplication experiment. An in- falling 
observer sends his information to the out- going direction (A); his information is also 
contained in the Hawking radiation (B). If the region beyond the inner horizon is 
regular, a second observer can see both A and B violating the no cloning theorem. 

black hole can see the information of the in-falling matter in the Hawking radiation. 
However, since the free-falling information is not affected by the Hawking radiation, 
two copies of the information exist, which violates the no cloning theorem of quantum 
mechanics. Here, black hole complementarity argues that this is not a problem since 
the two copies cannot be observed by a single observer. To check whether one observer 
can see both copies of the information or not, one may consider a gedanken experiment. 
One observer free-falls into the black hole and sends his information A in the out-going 
direction right after he crosses the event horizon. After the information retention time, 
a second observer sees information B in the Hawking radiation which includes A, and 
then goes into the black hole to get the information A. If this is possible, the second 
observer sees both A and B violating the no cloning theorem. However, in a neutral 
black hole, to deliver the information to the second observer before it collapses to the 
singularity, the first observer should send the signal with energy AE ~ exp M^, which 
is greater than the black hole mass itself, and hence this process is impossible. This 
argument also holds for charged black holes with the causal structure in Figure 17. 

However, if we assume sufficiently large A^, so that we can extend the causal 
structure of Figure 14 beyond the inner horizon, the duplication experiment will be 
possible, see Figure 18. Since the first observer can send a signal A in the out-going 
direction in time At ~ M, the signal requires energy AE ~ 1/M, which is sufficiently 
smaller than the black hole mass M, and hence this process is possible. Thus, the 
second observer can see both information A and B in the region beyond the inner 
horizon, violating black hole complementarity. For further discussions, see [36]. 
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6. Conclusion 

We constructed a spherically symmetric charged black hole model with a complex scalar 
field and a gauge field, and used the 1+1 dimensional approximation for the renormalized 
stress tensor to include Hawking radiation. Solving the equations numerically, we 
determined the Penrose diagram for dynamical charged black holes (Figure 14). We 
found a space-like inner horizon separated from the Cauchy horizon, which is parallel to 
the out-going null direction. Also, we demonstrated the transition from a charged black 
hole to a neutral black hole, and observed the generation of a space-like singularity and 
a second Cauchy horizon. 

Inner horizons have often been regarded as Cauchy horizons and null curvature 
singularities. However, we demonstrated that, in our model, the inner horizon is 
distinct from the Cauchy horizon, and it is not a curvature singularity from the general 
relativistic point of view; it is regular and penetrable. However, if the number of scalar 
fields are less than exp M^, the curvature functions blow up to greater than the Planckian 
cutoff beyond the inner horizon, so that we cannot extend physics beyond the inner 
horizon. Hence it is fair to say that a charged black hole has a "physical" space- 
like singularity even in dynamical cases (Figure 17). However, if we can assume an 
exponentially large number of massless fields, we may trust our results beyond the inner 
horizon. Then we found that strong cosmic censorship and black hole complementarity 
can be violated. Of course, the consistency of the assumption of an exponentially large 
number of massless fields should be checked within string theory. 
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Appendix A. Equations and integration schemes 

We use the gauge-invariant Lagrangian with a complex massless scalar field 4> and an 
electromagnetic gauge field ^4^ [2] : 



where Fab = ^b;a — ^a;b and e is the gauge coupling. From this Lagrangian we can derive 
the equations of motion for the scalar field and the electromagnetic field: 



>C = - (0;a + ieAa(l>) g"' (0., - ieAbcf>) - ^F^^F' 



(A.1) 




+ ieA" {2(j),a + ieAa(l)) + ieAa-bg (j) = 0, 

.6 - ie0 - ieAa^) + ie0 (0;„ + ieA^cj)) = 0. 



(A.2) 
(A.3) 
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Also, the energy-momentum tensor becomes 

Tab = \ {<P;a'<P;b + 0;a0;fe) 

1 — 1 

+ —FacFb' + e^AaAb(j)(j) + -Cgab. (A.4) 

We use the double-null coordinates [23], 

ds^ — — a'^{u, v) du dv -\- r'^{u, v) dO,^ , (A. 5) 

assuming spherical symmetry. Because of spherical symmetry and gauge symmetry, we 
can choose the gauge field as A/j, — (a, 0, 0, 0) with a single function a{u, v) [21]. 

Now, we will describe our numerical setup. We follow the notation of [19] [20] [21]: 
the metric function a, the radial function r, the Maxwell field a, and the complex 
massless scalar field s = \/47r0, and define 

h=—, d=—, / = r„, g = ry, w = Su, z = s„. (A.6) 

a a 

Then the Einstein and energy-momentum tensor components are 

Guu^ --{fu-Vh), (A.7) 
r 

Guv^^{^rU + a'' + ^fg), (A.8) 

G.. = --{g,-2gd), (A.9) 
r 

Gee = -^^(^du+^y (A. 10) 

Tuu — ^ [ww -|- iea{ws — ws) + e^a^ss] , i-^Al) 

= ^^z, (A. 13) 

Tee = / n {wz + zw) + ieaCzs - zs) + 
The scalar field and Maxwell field equations become 

0?C[ 

"fZu + jz + gw -|- iearz -|- ieags -|- ies— — = 0, (A. 15) 

^) ^'^[z-s-sz) =0, (A.16) 



(A.14) 



where 



?= — T- A.17 



can be interpreted as the electric charge. 
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We use the semi-classical Einstein equation, 

G^. = 87r (r° + (f^^j) (A.18) 

to include Hawking radiation. Spherical symmetry makes it is reasonable to use the 
1 + 1-dimensional results for (T^^)[24] divided by Anr'^ [20]: 

(Tl) -£rAhu- h') , (A.19) 

(T^) = (tn) = - (A.20) 

(t".> -^r^id.- d') , (A.21) 

with P = A^Z|i/127r, where N is the number of massless scalar fields and Ipi is the Planck 
length. 

Finally, we can list the simulation equations: 
(i) Einstein equations: 

- —{sz-^z)\, 
p 

= 2dg - rzz {dy - d'^) , (A.23) 

r 

fu — '^fh — rww — iear{ws — ws) — e'^a'^rss 

-^{h^-h'), (A.24) 

fg a'^ a^q^ P 
fv^Qu^ -p + -r-s d^. A.25) 



(ii) Maxwell equations: 

2^' 
e 

"2 



a. = ^, (A.26) 



2 



= - —-{sz-sz). (A. 27) 



(iii) Scalar field equations: 



fz gw iearz ieags ie 2 



= = T^OL qs. (A.28) 

Also, we have the definitions in Equation (A. 6). If we substitute Equation (A. 22) into 
Equation (A.25), then all equations contain only one derivative, except Equation (A.23) 
and (A.24). 

We can choose two different integration schemes corresponding to a choice of the 
equation. First, we can get r from g using Equation (A.23), a from d, and s from z. 
Second, we can get r from / using Equation (A.24), a from /i, and s from w. We mainly 
used the first integration scheme. If we call the radial function of the first scheme r(t,) 
and the radial function of the second scheme r(„), then r(i,) and r(„) should be the same. 
We compared them to check the consistency of our simulation. 
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We solved these equations using the second order Runge-Kutta method [25] 
following [21]. Equation (A. 22) has a singularity at r = \/P and we may regard it 
as the central singularity due to the semi-classical approximation. For consistency y/P 
should be sufficiently smaller than the size of the black hole in the simulation. 

Appendix B. Initial conditions 

We need initial conditions for all functions {a, r, g, f, h, d, s, w, z, a, q) on the initial 
u — Ui and v — Vi surfaces, where we set Ui = vi = 0. 

We have gauge freedom to choose the initial r function. Although all constant u 
and V lines are null, there remains freedom to choose the distances between these null 
lines. Here, we choose r(0,0) = tq, f{u,0) = Vuo, and g{0,v) = r^o, where r„o < 
and r„o > so that the radial function for an in-going observer decreases and for an 
out-going observer increases. 

We use a shell-shaped scalar field, and hence its inside is not affected by the 
shell. Thus, we can simply choose q{u, 0) = a{u, 0) = and a{u, 0) = 1. Also, 
s{u, 0) = w{u, 0) = h{u, 0) = holds. Then, since the mass function, defined in 
Equation (8), should vanish at w = v = 0, it is convenient to choose r^o = ~l/2 
and r„o = 1/2. 

We need more information to determine d, g, and z on the v = surface. We get d 
from Equation (A. 22), g from Equation (A. 25), and z from Equation (A. 28). 

We can choose an arbitrary function for s(0, f ). For example, to make a black hole, 
we use 

s(0, v) = A sin^ (^^'^^ (^^^^i") ^^'^"^ 

for < V < V{ and otherwise s(0, v) — 0, where Vf is the width of the pulse and A is 
the amplitude. Then we obtain z{0,v). Note that, according to Equation (A. 27), if we 
represent the scalar field as s = |s| exp(iQ), then — er^|s|^Q^ holds, and the sign of 
the exponent determines the sign of charge. 

Also from Equation (A. 23), we can use d = rz'z/2g on the u = surface, since we 
can assume that there is no Hawking effect on the initial surface, and we get d{0, v). By 
integrating d along v, we get a{0,v). 

We need more information for h, /, w, a, and q on the u — surface. We obtain q 
from Equation (A. 27) and a from Equation (A. 26). Then, we get h from Equation 
(A.22), / from Equation (A.25), and w from Equation (A.28). This finishes the 
assignments of initial conditions. 

We choose Tq = 10 and V{ — 20, leaving the three free parameters (e, P, A). 

Appendix C. Convergence and consistency checks 

To check the convergence of our simulations, we compared the results for (e = 0.1, P — 
0.1, A — 0.25) using different step sizes: 1, 2, and 4 times finer. In Figure CI, we see 
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0.1 r u=15 I ; u=10 



0.001 r / /'^ -- 

0.0001 r /■'/'/ - 

le-05 1 1 1 1 1 

10 20 30 40 50 60 10 20 30 40 50 60 

V V 

Figure CI. Plots of errors with different step sizes. Here, we plot |r(i) — ?'(2)|/''(2) 
(solid curves) and 4|r(2) — ?'(4)|/?'(4) (dashed curves) along a few constant u lines, 
where is calculated in an n times finer simulation than r(^iy This shows that our 
simulation converges to second order. 

Figure C2. |r(^,) — ?'(u) along a few constant u lines, where r(y) is calculated 

by integrating g using Equation (A. 23), and r(-„-) is calculated by integrating / using 
Equation (A.24). 




that the difference between 1 and 2 times finer is 4 times the difference between 2 and 
4 times finer, and thus our simulation converges to second order with errors < 0.1%. 

To check the consistency, we used two independent integration schemes as 
mentioned in Appendix A. Figure C2 shows that the difference between the two schemes 
is < 1%. 

Near the outer horizon, the radial function r changes rapidly causing the slow 
convergence shown by the spikes in Figure CI and C2. To avoid this problem, we chose 
the step size Au adaptively so that the ratio ruAu/r was constant [21]. 

Finally, we checked whether our simulation gives the physically correct picture for 
a neutral black hole [20]. Figure C3 and C4 are contour diagrams for r and the outer 
horizon r„ = for (e = 0, P = 0.1, A = 0.25). As the matter collapses, the outer horizon 
bends in a space-like direction. After the matter collapse ends, the outer horizon bends 
in a time-like direction due to the Hawking radiation. These phenomena are consistent 
with the properties of local horizons [29] . Moreover, one can see the space-like central 
singularity at r = y/P. This confirms that our simulation gives the correct result for a 
neutral black hole (Figure C5) [37]. 
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